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ABSTRACT 


Using  the  geometrical  optics  (GO)  and  physical  optics  (PO)  approximations, 
a  correct,  complete,  ready-to-use  formula  is  derived  for  the  backscatter  (monostatic) 
polarization  scattering  matrix  (PSM)  of  the  perfectly  conducting  right  dihedral  at 
arbitrary  incidence  angle.  The  absence  of  such  a  result  from  the  literature  is  sur¬ 
prising  given  that  the  dihedral’s  PSM  is  needed  in  many  applications,  such  as  in  the 
calibration  of  polarimetric  radars,  including  synthetic  aperture  radars  (SAR),  in  the 
generation  of  simulated  polarimetric  radar  imagery,  and  in  automatic  target  recog¬ 
nition  (ATR).  Because  the  new  results  provided  are  important  to  many  researchers 
who  may  not  be  experts  in  electromagnetic  theory  (as  is  often  the  case  for  the 
computer- vision  researchers  working  on  ATR),  the  report  is  relatively  self-contained 
and  takes  the  reader  from  the  definitions  of  PSMs  and  complex  radar  cross-sections, 
through  the  mathematical  formulation  of  Huygen’s  Principle,  the  combined  use  of 
GO  and  PO,  and  changes  of  polarization  bases,  to  the  derivation,  discussion,  and 
simplification  of  the  dihedral’s  PSM. 
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1.  INTRODUCTION 


The  right-angle  dihedral  (DI)  corner  reflector  is  an  important  reference  device  for  calibrating 
radar  systems.  It  is  often  used  in  conjunction  with  the  right-angle  trihedral  (TRI)  corner  reflector, 
the  scattering  properties  of  which  are,  to  some  degree,  complementary. 

1.  The  strongest  returns  from  a  DI  are  due  to  double-bounce  interactions,  in  which  the 
incident  wave  successively  hits  both  faces  of  the  DI  before  being  returned  towards  the 
source  (thus,  DIs  are  known  as  double-bounce  or  even-bounce  scatterers).  In  contrast, 
the  strongest  returns  from  a  TRI  are  due  to  triple-bounce  interactions,  in  which  the 
incident  wave  hits  all  three  faces  of  the  TRI  in  succession  (thus,  TRIs  are  known  as 
triple-bounce  or  odd-bounce  scatterers). 

2.  The  DI  produces  a  strong  return  for  simultaneous  transmit  right  (or  left)  circular 
polarization  and  receive  right  (or  left)  circular  polarization.  In  contrast,  the  TRI 
produces  a  strong  return  for  opposite  transmit  and  receive  circular  polarizations. 

3.  The  backscatter  pattern  of  the  DI  is  broad  in  a  plane  perpendicular  to  its  crease  and 
narrow  in  planes  containing  it.  In  contrast,  the  backscatter  pattern  of  the  TRI  is 
relatively  narrow. 

Even  though  these  properties,  as  well  as  others,  are  well  documented,  the  literature  does  not 
readily  provide  formulas  for  the  backscatter  (or  monostatic)  polarization  scattering  matrix  (PSM) 
of  elementary  scatterers  (such  as  the  dihedral  and  trihedral)  at  arbitrary  incidence  angles. 

The  absence  of  such  formulas  may  result  from  the  fact  that  fully  polarimetric  radars  are 
relatively  new.  The  recent  availabiltity  of  fully  polarimetric  SAR  imagery  (e.g.,  from  the  ERIM  [1], 
JPL  [2],  and  MIT  [3]  airborne  sensors)  has  led  to  a  resurgence  of  interest  in  polarimetry  and  to  a 
quest  for  analytical  formulas  describing  the  monostatic  PSMs  of  elementary  scatterers  at  arbitrary 
incidence  angles. 

Accurate  PSM  formulas  for  elementary  scatterers  are  needed  by  a  class  of  computer  codes  pre¬ 
dicting  the  appearance  (or  signature)  of  complex  targets  in  polarimetric  range  profiles  or  range/cross¬ 
range  images.  For  example,  SarTool  [4]  first  decomposes  CAD  models  of  complex  targets  into  their 
constituent  elementary  scatterers  and  then  essentially  add  their  properly  transformed  PSMs  to 
produce  the  PSM  of  the  targets.  In  the  case  of  range  profiles  and  range/cross-range  images,  the 
calculation  is  repeated  for  each  point  in  the  profile  or  image. 

Simple,  closed-form  expressions  for  the  PSMs  of  elementary  scatterers  are  also  needed  by 
computer- vision /image-understanding  researchers  developing  automatic  target  recognition  (ATR) 
algorithms.  In  one  line  of  research,  bright,  contrasted,  point-like  features  are  extracted  and  labeled, 
e.g.,  as  being  “dihedral-like” .  Generally,  these  labels  are  found  using  techniques  based  on  the  early 
work  of  Huynen  [5]  (e.g.,  see  Cameron  and  Leung  [6]).  However,  much  work  remains  to  be  done  to 
solve  this  difficult  “inverse  problem.” 
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A  number  of  techniques  can  be  used  to  predict  the  PSMs  of  elementary  scatterers,  but  com¬ 
bining  geometrical  optics  (GO)  [7,  p.  869]  and  physical  optics  (PO)  [7,  p.  870]  provides  the  simplest 
possible  formulas  still  having  a  reasonable  degree  of  accuracy.  In  the  case  of  the  dihedral  and  tri¬ 
hedral,  the  author  is  not  aware  of  any  report  providing,  by  one  method  or  another,  ready-to-use, 
accurate,  closed-form  PSM  formulas  at  arbitrary  incidence  angle.  Blejer  [8,  9]  comes  close  to  pro¬ 
viding  such  formulas,  but  the  results  he  gives  should  be  questioned  because  of  likely  errors  and 
omissions  in  the  paper  by  Corona  et  al.  [10]  from  which  these  results  are  derived. 

The  first  goal  of  the  present  report  is  to  correct  the  above  problems  for  the  dihedral,  thereby 
providing  an  accurate  set  of  formulas  for  the  dihedral’s  monostatic  PSM  resulting  from  combined 
GO  and  PO.  The  second  goal  is  to  make  the  report  accessible  to  researchers  who  need  such  formulas 
and  are  not  experts  in  electromagnetic  (EM)  theory,  as  is  often  the  case  in  the  ATR  and  computer- 
vision  communities. 

Now,  a  brief  chronological  review  of  relevant,  prior  work  on  the  analysis  of  the  dihedral’s  PSM 
is  provided. 

Knott  [11]  considers  an  obtuse- angle  dihedral,  with  illumination  perpendicular  to  the  crease 
and  with  arbitrary  polarization.  Single-bounce  (SB)  and  double-bounce  (DB)  returns  are  handled 
by  using  PO  to  deal  with  scattering  by  the  last  surface  encountered  (in  the  SB  and  DB  cases)  and 
GO  to  deal  with  reflections  from  the  first  surface  in  the  DB  case.  The  goal  of  the  paper  is  to  study 
the  effect  of  the  dihedral  angle  on  radar  cross  section. 

Corona  et  al.  [10]  extend  Knott’s  work  to  arbitrary  incidence.  Even  though  their  results  are 
correct  in  the  case  of  normal  incidence,  they  are  believed  to  be  erroneous  and  incomplete  in  the 
general  case. 

Blejer  [8,  9]  used  Corona’s  incorrect  and  incomplete  results  to  derive  the  PSMs  of  the  dihedral 
and  trihedral.  These  formulas  are  apparently  used  in  the  prediction  code  SarTool  [4]. 

Griesser  et  al.  [12]  consider  a  problem  similar  to  Knott’s  (incidence  normal  to  crease;  vertical 
and  horizontal  polarizations),  but  they  deal  with  arbitrary  dihedral  angles,  as  well  as  with  single-, 
double-,  and  triple-bounces,  single  refractions,  and  reflection- diffractions,  using  various  combina¬ 
tions  of  GO,  PO,  and  PTD  (physical  theory  of  diffraction). 

Corona  et  al.  [13]  use  a  generalization  of  PO  to  study  the  effect  of  surface  loading  in  the  case 
of  orthogonal  incidence. 

Anderson  [14]  considers  vertical  and  horizontal  polarizations  at  arbitrary  incidence  angles 
on  dihedrals  with  angles  no  smaller  than  60  deg,  and  uses  PO  to  compute  single-,  double-,  and 
triple-bounce  terms.  However,  the  integrated  forms  of  key  surface  integrals  are  not  given  in  the 
paper. 

Atkins  and  Shin  [15]  propose  a  new  approach  for  dealing  with  multiple  bounces  between  two 
possibly-disjoint  polygonal  plates  (a  special  case  of  which  is  the  dihedral),  where  they  avoid  using 
GO  for  dealing  with  intermediate  reflections. 
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This  report  explores  in  detail  the  derivation  of  the  PSM  of  the  (unloaded)  right  dihedral 
under  the  approximations  of  GO  and  PO.  As  already  indicated,  it  is  written  to  make  the  informa¬ 
tion  accessible  to  the  nonexpert,  in  particular  to  researchers  in  computer- vision  and  ATR.  Thus, 
this  report  begins  with  a  brief  review  of  fundamental  concepts,  such  as  PSMs,  complex  radar 
cross  sections,  changes  of  polarization  basis,  GO,  and  PO.  Then,  a  general  expression  for  the  field 
backscattered  by  the  dihedral  is  derived,  avoiding  the  premature  combination  of  SB  and  DB  terms. 
An  important  dichotomy  arising  in  the  integration  of  a  DB  surface  integral  is  carefully  considered. 
Finally,  the  PSM  of  the  dihedral  is  derived,  the  relative  contributions  of  the  SB  and  DB  terms 
are  compared,  and  simplified  formulas  are  derived,  which  agree  with  the  expression  used  almost 
exclusively  by  those  working  the  inverse  problem  (e.g.,  Huynen  [5]  and  Cameron  and  Leung  [6]). 
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2.  BACKGROUND 


2.1  Fields 

Consider  time-harmonic,  i.e.,  monochromatic,  EM  fields  with  angular  frequency  u>  =  2irc/\, 
where  c  is  the  speed  of  fight  and  A  the  wavelength  of  interest.  Using  a  time-dependence  e3ujt ,  the 
electric  field  at  any  point  fp  can  be  written  as  [16,  p.  46] 

E(rp,t)  =  ne{E(rP)  (1) 

This  report  focuses  on  the  complex,  time-independent  field  E(fp),  also  known  as  a  phasor.  The 
physical  field  can  be  retrieved  at  any  time  through  Equation  (1).  These  remarks  also  apply  to  the 
magnetic  field  H. 

In  the  classical  treatment  of  scattering,  the  wave  incident  on  the  scatterer  is  assumed  to  be 
a  plane  wave  and  the  scattered  wave  is  an  outgoing  spherical  wave  that,  in  the  far  zone  of  the 
scatterer,  can  be  treated  locally  as  a  plane  wave. 

Plane  wave  are  characterized  by  their  propagation,  or  wave,  vector  k  =  kk,  where  k  —  2tt /A 
and,  in  general,  x  denotes  the  unit  vector  corresponding  to  x.  In  the  case  of  linear  polarization 
characterized  by  a  polarization  p  orthogonal  to  the  propagation  direction  k,  the  complex  field 
E(rP)  =  Ep(rP)  is  [16,  p.  54] 


Ep(fP)  =  p  (2) 

where  Ep  is  the  complex  magnitude  of  the  field  at  fp  =  0.  The  case  of  arbitrary  polarization  is 
considered  below. 

2.2  Scattering 

Consider  the  scattering  geometry  of  Figure  1.  The  scatterer  is  assumed  to  be  located  near 
the  origin  0  of  the  (a:,  y ,  z)  coordinate  system,  and  the  sensor  is  assumed  to  be  located  at  a  point  S 
with  spherical  coordinates  (r,  <f>).  With  each  S,  one  associates  a  unique  coordinate  system  defined 

by  the  unit  vectors  f,  8,  and  r  is  defined  by  the  vector  US  =  f;  $  is  perpendicular  to  the  plane 
defined  by  S  and  the  2- axis,  and  oriented  in  the  sense  of  increasing  (p:  and  6  is  perpendicular  to 
the  plane  ( f ,  <j>)  and  oriented  in  the  sense  of  increasing  6,  with  the  result  that  the  frame  (f ,  8 ,  4>)  is 
right-handed.  Of  course,  the  sensor  position  is  given  by  f  —  rf. 
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Figure  1.  Scattering  geometry. 


Consider  an  EM  field  transmitted  from  S  towards  0.  If  the  scatterer  is  in  the  far  zone  of 
the  transmitting  antenna,  this  field  can  he  treated  as  a  plane  wave  traveling  in  the  k  —  kt  =  -f 
direction.  In  the  case  of  linear  polarization  along  p  (perpendicular  to  f),  the  “transmitted”  complex 
field  E^rp)  can  be  written  as  in  Equation  (2).  In  the  case  of  arbitrary  polarization,  Et(rp )  can 
be  described  in  terms  of  its  components  in  a  polarization  basis  ( AB )  characterized  by  its  basis 
vectors  h a  and  hp.  A  particularly  convenient  basis  for  the  intended  analysis  is  the  (HV)  basis  (H 
=  horizontal,  V  =  vertical),  the  basis  vectors  of  which  are  the  vectors  <p  and  9  of  Figure  1,  i.e., 
hp  =  4>  and  hy  —  6.  Other  choices  of  basis  are  discussed  later.  In  the  (HV)  basis,  one  has 

Ei(fp)  =  Et(fP)  +  E|(fp)  (3) 

with  [see  Equation  (2)] 

=  ElP~ik"fF 
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(4) 


where  pt  is  either  4>  or  9 .  Equation  (3)  can  thus  be  written  as 

E*(rP)  =  E1  e~^tIp  (5) 

with 

Et  =  Et24>+EtJ.  (6) 

q>  u 

The  magnitudes  and  phases  of  the  complex  numbers  and  £|  fully  determine  the  polarization 

characteristics  of  the  transmitted  field  E*(rp).  Note  that  the  plane  wave  described  by  E\fP )  is  a 
local  approximation  to  a  spherical  wave  originating  from  S,  so  that  M  and  Ej  should  include  the 

l/r  attenuation  associated  with  the  distance  r  from  S  to  0.  However,  this  report  treats  E\fP )  as 
a  true  plane  wave  and,  thus,  E |  and  Ej  as  (complex)  constants. 

The  currents  induced  on  the  surface  of  the  scatterer  by  the  incoming  wave  produce  a  scattered 
EM  field.  In  the  far  zone  of  the  scatterer,  this  field  is  effectively  a  spherical  wave  originating  from 
0  that  can  be  approximated  locally  by  a  plane  wave.  Limiting  the  analysis  to  the  backscattering 
case  (i.e.,  to  transmit  and  receive  sensors  co-located  at  S),  the  “received”  E-field  can  be  expressed 
as 


Er(f)  =  Er-(f)P  +  El(f)9. 


(7) 


Note  the  use  of  f  (instead  of  rP)  since  one  is  only  interested  in  the  fields  at  the  receiver.  Also  note 
that  the  f-dependencies  have  not  been  factored  out  [in  contrast  to  what  was  done  in  Equations  (5) 
and  (6)  with  regard  to  fP]. 


2.3  Polarization  Scattering  Matrix 


It  is  well-known  [17,  18,  19]  that  the  complex  components  Er~(f)  and  Er-(r)  of  the  received 
field  Er(f)  are  linearly  related  to  the  complex  components  E|  and  Ej  characterizing  the  value  of 
the  transmitted  field  at  the  origin  0,  i.e.,  of  E\fP  =  0).  This  linear  relation  is  written  as 


( *F> 

v  m?) 


/ 

a(k, r ) 

\ 


(8) 
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or 


Er(r)  =  a(k ,  r)  S  El 


(9) 


where  the  use  of  the  scale  factor  a(k,  r )  reflects  the  lack  of  universal  agreement  on  the  definition 
of  the  (polarization)  scattering  matrix  5.  Typical  values  for  a(k,  r)  are 


1, 


e~ikr 

kr 


(10) 


Note  the  slight,  inconsequential  abuse  of  notation  for  El  and  ET(r).  Indeed,  these  symbols 
denote  3-D  vectors  [in  the  ( x,y,z )  frame]  in  Equations  (6)  and  (7),  and  2-D  vectors  [in  the  (<M) 
frame]  in  Equation  9. 

In  practice,  the  elements  Sprpf  of  S  (with  pT  and  pt  taking  either  of  the  values  <j>  or  9)  can 

be  found  as  follows.  Consider  a  transmitted  field  that  is  linearly  polarized  along  pt  =  p  or  pt  =  6 

and  has  some  arbitrary  complex  magnitude  Eq.  Such  a  field  can  be  described  by  Equation  (4)  with 

Ei  -  Eq.  From  Equation  (6),  it  follows  that 
1 Pt 

E^EoPt.  (n) 


Using  this  expression  in  conjunction  with  Equations  (8)  and  (9),  one  finds 


~  r)  Sprpt  Eo 


(12) 


so  that 


Sprpt  a.{k ,  r ) 


EpSf) 


Eq  ' 


(13) 


Because  Ep  is  the  component  along  pr  of  the  field  Eph)  received  when  the  transmitted  field  has 
complex  amplitude  Eq  and  polarization  pt,  one  can  be  rewrite  (for  future  use)  Equation  (13)  as 


1  Ept{f).pr 

a(k,r )  E0 


(14) 
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2.4  Complex  Radar  Cross-Sections 

Because  of  the  lack  of  agreement  on  the  definition  of  a(k,  r),  it  is  useful  to  first  find  the 
scattered  complex  radar  cross  sections  (C-RCS)  yf°prpt  for  tiie  various  combinations  of  pt  and  pr, 
and  then  express  each  S^pt  in  terms  of  the  corresponding  ^p~ft  (note  that  ^prpt  is  the  symbol 
denoting  each  C-RCS).  The  more  familiar  radar  cross  sections  (RCS)  cprpt  are  obtained  by  taking 
the  square  of  the  magnitude  of  y/afirpi  ■ 

The  agreed-upon  definition  for  the  C-RCS  is  [8] 


lim  2\/7rr 

r—voo 


Ept(r).pr 

E0 


eikr 


(15) 


In  all  cases  of  interest  here,  the  limit  in  Equation  (15)  can  be  ignored,  so  that  the  desired  relationship 
is 


SprPt  ~  ^  \l°PrPt 


(16) 


where 


,—ihr 


2y/tr  a(k,  r) 

The  values  of  /3(fc,r)  for  the  three  common  choices  in  Equation  (10)  are 


(17) 


ikr 


2yph rr  5  2 ■spit  ’  2a/7t 


(18) 


Equation  (15)  can  be  used  to  define  a  C-RCS  matrix  yfa  analogous  to  the  scattering  matrix 
S  in  Equations  (8)  and  (9),  i.e., 


/ 


\fo 


77  \ 


fu  \r&) 


(19) 
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Clearly 


5  =  /3(k,  r)y/a. 


(20) 


2.5  Change  of  Polarization  Basis 

The  previously-selected  basis  vectors  hR  =  cf>  and  hy  —  0  of  the  linear  basis  (HV)  are  true 
3-D  space  vectors  with  real-valued  components.  However,  in  general,  the  basis  vectors  hA  and  hB 
defining  a  basis  ( AB )  may  be  three-component  vectors  with  complex- valued  components.  A  typical 
example  is  the  left-right  circular  basis  (LR)  with  basis  vectors  hi  and  Kr  related  to  Hr  and  hy  by 


(  hL  \  =  _L  (  1  j 

\  hR  )  ~  V2  \  1  -j 


or,  symbolically, 

h(LR)  =  T(HV-+LR)h(HV) 


where 


1 

R«v-W  =  yi 


1 

1 


j 


(21) 


(22) 


(23) 


and,  in  general,  h(AB )  is  a  “vector,”  the  elements  of  which  are  the  unit  vectors  hA  and  hB  associated 
with  the  basis  (AB).  In  the  general  case  of  two  bases  (AB)  and  (A'B'),  one  writes 


h(A'B’)  -  T(AB^A'B')h(AB)- 


(24) 


Although  the  transformation  T(jjv-+lR)  ln  Equation  (23)  appears  to  be  the  one  most  often 

used  to  relate  (HV)  and  (LR),  this  transformation  is  not  universally  agreed  upon  [20,  p.  17],  [8, 
p.  4]. 
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The  effect  of  a  “change  of  basis”  on  the  scattering  matrix  is  a  delicate  subject  treated  in 
detail  in  Appendix  A,  where  it  is  shown  that  a  dual  change  of  basis  is  actually  used,  with  one 
change  for  the  transmitted  field,  and  one  separate,  related  change  for  the  received  field.  If  the 
scattering  matrix  is  given  by  S^ab)  when  the  transmitted  field  E*  is  expressed  in  basis  (AS),  then 
Equation  (A.l)  and  (A.2)  show  that  the  scattering  matrix  S(A'B‘)  corresponding  to  E*  expressed 
in  basis  ( A'B ')  is 

S(a’B')  =  UTS(ab)U  (25) 


with 


U  =  Tt  (26) 

where,  in  general,  AT  denotes  the  transpose  of  A.  Obviously,  Equation  (25)  can  also  be  written  as 
S(a'B')  =  TS(ab)Tt.  (27) 

Carefully  note  that  one  has  not  said  in  which  basis  the  received  field  ET  is  expressed  and  that 
Equation  (25)  cannot  generally  be  derived  from  Equation  (8)  by  performing  a  true  change  of  basis 
(these  points  are  discussed  in  Appendix  A). 

For  calculation  purposes,  observe  that  U  and  T  are  both  unitary  matrices  [see  Equations  (A.9) 
and  (A. 10)].  Finally,  for  future  reference,  note  that  the  relation  between  the  linear-basis  scattering 

matrix 


S(HV)  - 


Shh  Shv 
Shv  Svv 


and  the  corresponding  circular-basis  scattering  matrix 


(28) 


S(LR) 


Sll  Slr 
Slr  Srr 


(29) 
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is  given  by 


l/l  3 


J(^)  =  2  l  i  -j 


Shh  Shv 
Shv  Svv 


1  1 

3  -3 


(30) 


or 


S(lR)  - 


Shh-Svjl  +  j  sHV  sw+Sxk 

Shh-Svv  _jsHV  j 


Srw+Svv 

2 


(31) 


Clearly,  a  similar  formula  can  be  written  for  the  C-RCS  matrix  y/a. 

2.6  Huy  gen’s  Principle 

The  key  formula  for  computing  the  field  Ept  appearing  in  Equations  (14)  and  (15)  is  the 
mathematical  expression  of  Huygen’s  principle  [16,  p.  399,  hut  with  e*"1],  i.e., 


(32) 


where  a  is  the  permeability  of  the  medium  (assumed  isotropic),  S'  the  surface  of  the  scatterer  (with 
points  r').  and  JsPt(?)  the  surface  current  induced  on  S'  by  a  transmitted  wave  wrth  polarisation 

Pt- 


The  dyadic  operator  [I  -  ff ]  applied  to  some  arbitrary  vector  a  should  be  interpreted  as 


[I  -  ff]a  -a-  r(r.a ) 


(33) 


which  has  the  effect  of  removing  the  component  of  a  along  r,  thereby  yielding  the  projection  of  a 
in  the  plane  ( 0,<j> ),  i.e., 


[j  _  fr]a  =  4>(4>-ai)  +  0(9.a). 


(34) 
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The  dyadic  operator  can  thus  be  rewritten  as 


J  —  r?  =  <p(j>  +  69.  (35) 

Note  that  e~ikr [I  -  fr]  is  a  mathematical  expression  of  the  fact  that  Ept(r] )  is  an  outgoing 
spherical  wave,  which  is  precisely  the  assumption  made  previously  in  the  discussion  leading  to 
the  definition  of  the  scattering  matrix  S.  In  fact,  Equation  (32)  can  be  used  to  obtain  explicit 
expressions  for  the  4>  and  9  components  of  Epjr)  that  appear  in  Equation  (7),  but  this  is  not 
further  explored  here. 

2.7  Geometrical  and  Physical  Optics  Approximations 

One  needs  to  relate  the  surface  current  J s^pt{rP)  to  the  transmitted,  linearly  polarized  (along 
pt )  wave  characterized  by  [see  Equation  (4)] 

Elpp)  =  E‘pt  p,  =  E<pt  p„  (36) 

For  compatibility  with  Equations  (14)  and  (15),  is  set  to  E0  [see  Equation  (11)],  so  that 
Equation  (36)  becomes 

£!,( fF)  =  Eo  e-ihfp  P,  =  E„  eikf  fp  p„  (37) 

To  find  jsptir'),  one  also  needs  the  magnetic  field  f?|t  associated  with  [16,  p.  55],  i.e., 

HUip)  =  —  k,  x  EUfF)  =  -EEo f  x  pt-  (38) 

Pt  up  pt  up 

The  counterpart  of  Equation  (37)  is  thus 


HUrp) 


pt' 


H0e~fk^Pat  =  Hoeikf-rPat 


(39) 
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with  the  complex  value  Hq  given  by 


Ho 


OJfl 


and 


at 


6  if  pt  =  <i> 
~4>  if  Pt  =  <?• 


(40) 


(41) 


For  a  number  of  elementary  scatterers  of  interest  (e.g.,  the  dihedral  considered  later  in  the 
report),  it  is  necessary  to  consider  not  only  single-bounce  interactions,  but  also  multi-bounce  inter¬ 
actions  (e.g.,  double-bounce  in  the  case  of  the  dihedral).  Here,  one  follows  the  customary  practice 
of  using  PO  to  compute  the  field  backscattered  by  the  last  surface  participating  in  a  particular 
interaction,  and  of  using  GO  to  handle  scattering  by  intermediate  surfaces.  The  following  analysis 
relies  on  the  assumption  that  all  scattering  surfaces  are  perfect  conductors. 

First,  the  magnetic  field  incident  on  the  last  scattering  surface  (S')  is  derived  from 

Hpt{fp)  and  the  scatterers  geometry  through  the  laws  of  GO.  Second,  one  uses_the  PO  approxima¬ 
tion  to  relate  the  surface  current  JSpt(r')  appearing  in  Equation  (32)  to  if|t(r'),  i.e.,  [7,  Equation 
(38)], 


Js,pt(r')  ~ 


2  n(r')  X  Hpt(r')  on  illuminated  part  of  S' 
0  on  shadowed  part  of  S' 


where  n(r')  is  the  unit 


normal  to  S'  at  r'. 


(42) 
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3.  BACKSCATTERED  E-FIELD  FOR  DIHEDRAL 


In  this  section,  compact,  nonintegral  expressions  are  derived  for  the  field  Ept(r)  [given  by 
Equation  (32)]  backscattered  from  a  perfectly  conducting  dihedral  sized  and  positioned  as  indicated 
in  Figure  2.  Recall  that  the  subscript  pt  indicates  that  the  field  corresponds  to  a  transmitted  plane 
wave  linearly  polarized  along  pt.  This  wave  is  characterized  by  the  E-field  of  Equation  (37)  or, 
equivalently,  by  the  H-field  of  Equation  (39). 

209852-1 


Figure  2.  Dihedral  geometry. 
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3.1  Single-Bounce  Integral  Expression 

The  dihedral  exhibits  two  distinct  single-bounce  (SB)  interactions.  The  first  corresponds  to 
backscattering  of  the  transmitted  wave  by  the  surface  Si  of  the  side  normal  to  y,  and  the  second  to 
backscattering  of  the  same  wave  by  the  surface  S2  of  the  side  normal  to  x.  Both  interactions  can 
be  treated  jointly  by  denoting  the  backscattering  surface  by  S  and  its  normal  towards  the  sensor 
by  h. 

In  the  SB  case,  the  field  appearing  in  Equation  (42)  is  simply  given  by  the  transmitted 
H-field  Hpt,  i.e.  [see  Equation  (39)], 


HU?)  =  Ho  e~iktT  at  =  H0  elKT:r  at. 


ikr.r 


pt 


Using  Equations  (32)  and  (42),  one  finds  that  the  SB  contribution  due  to  S  alone  is 


Ef, t(f)  =  ~tkEllr [/  "  X  ht)  Js 


e2ikr.T'  dS 


(43) 


(44) 


where  r'  denotes  points  in  S. 

3.2  Double-Bounce  Integral  Expression 

The  dihedral  exhibits  two  distinct  double-bounce  (DB)  interactions.  The  first  corresponds  to 
backscattering  by  S2  of  the  portion  of  the  transmitted  wave  that  is  reflected  by  Si  and  intercepted  by 
52,  and  the  second  is  obtained  by  reversing  the  roles  of  Si  and  S2-  Both  interactions  can  be  treated 
jointly  by  denoting  the  reflecting  surface  (the  first)  and  its  normal  by  SA  and  hA,  respectively,  and 
the  backscattering  surface  (the  second)  and  its  normal  by  SB  and  nB,  respectively.  The  part  of  Sb 
illuminated  by  reflections  from  SA  is  denoted  by  SAb- 

In  the  DB  case,  the  field  j?|t  appearing  in  Equation  (42)  is  obtained  by  reflecting  the  trans¬ 
mitted  H-field  H\t  [given  by  Equation  (39)]  essentially  according  to  GO: 

1.  The  wavevector  h  is  obtained  from  kt  by  requiring  that  the  reflection  on  SA  reverse 
the  sign  of  the  component  of  kt  corresponding  to  hA  (necessarily  either  £  or  y),  i.e., 

h  =  kt~  2(kt.nA)nA.  (45) 

2.  The  polarization  vector  d;  is  obtained  from  at  by  imposing  the  boundary  conditions 
for  the  H-field  on  the  perfectly  conducting  surface  SA,  be.,  that  (a)  the  tangential 
components  of  the  incident  and  reflected  H-fields  be  equal,  and  that  (b)  the  sum  of 
their  normal  components  be  zero.  One  finds 

a,i  =  at  -  2{at.hA)nA.  (46) 
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Therefore 


Jj-Jr')  =  )f‘ ■<!''[«:  -  2(a,.nA)nA] 

=  Boe* ^)*Al-r_'|a(  _  2(a«.ftx)nA].  (47) 

Again,  using  Equations  (32)  and  (42),  one  finds  that  the  DB  contribution  due  to  surfaces  SA 
and  Sb  encountered  in  that  order  is 

£?aSs( f)  =  -ikEf~'tri/-ff](nBxlat-2(d,.nA)nA])  f  iSAB(*$) 

Pt  Zitr  JSAb 

where  r'  denotes  points  in  SAB- 

3.3  Complete  Integral  Expression 

The  Equations  (44)  and  (48)  for  the  SB  and  DB  contributions  can  be  combined  to  get  the  total 
E-field  backscattered  by  the  dihedral,  i.e.,  [notice  the  similarity  with  Equation  (5)  in  Anderson’s 
paper  [14]], 

)  (49) 

with 


Ps  =  ( r.h )  [  e2ikf'?dS 

J  s 

(50) 

-S  [I-rf](hxat) 

HPt  ~  r.h 

(51) 

PSaSB  =  (r.hB)  [  e^[r-{f.nA)hA].P dsAB 

JsAB 

(52) 

qSaSb  [I  ~  rr](hB  x  [df  -  2 (at.hA)hA}) 

^Pt  r.nB 

(53) 

The  reason  for  introducing  the  apparently  superfluous  factors  ( r.h )  and  ( r.hB )  will  become  clear 
later. 


Ep,(f)  =  - 


ikEoe 


—ikr 


2  irr 


PSlQSpl  +  PS2QSp2t  +  PSlS2QSpf2  +  PS2‘SiQpl 
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In  the  next  sections,  one  derives  compact,  closed-form  expressions  for  the  four  quantities 
above.  To  do  so,  one  will  express  r  in  terms  of  its  components  lx,  ly,  and  lz ,  which  are  also  the 
direction  cosines  of  the  sensor  location,  i.e., 


T  —  (lx,  ly,  lz) 

with 

lx  —  sin#  cos<?!> 

<  ly  =  sin  6  sin  </> 
lz  =  cos  6. 

3.4  Simple  Expressions  for  P^1  and  P^2 

When  S  =  Si,  one  has  h-y  and  r'  =  ( x,0,z ),  so  that  Equation  (50)  becomes 


i.e., 

pS 1  =  ably  sine  k'alx  sine  k'blz 


where 


sm7ra: 
smex  =  - 

7 TX 


and 


(54) 


(55) 


(56) 


(57) 


(58) 


(59) 


When  S  =  S2,  one  has  h  =  x  and  r'  =  (0 ,y,z),  so  that  Equation  (50)  similarly  becomes 


pS 2  =  ablx  etka^y  sine  k'aly  sine  k'blz. 


(60) 
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3.5  Simple  Expressions  for  Q j?1  and 

One  begins  by  rewriting  the  numerator  A  of  the  Equation  (51)  for  Qpt  as 

A  =  [/  -  ff](n  x  at)  (61) 

=  8(6. (h  x  at ))  +  4>(4>-(h  x  at ))  (62) 

=  ^(n.(ai  x  0))  +  4>(h.(at  x  p))  (63) 

where  one  has  used  Equation  (34)  and  the  cyclic  property  of  the  mixed  product,  i.e.,  a.(b  X  c)  = 
b.(c  x  a). 

1.  When  pt  =  <£,  Equation  (41)  gives  dj  =  #.  Then  Equation  (61)  reduces  to  (h.r)<j>, 
and  Equation  (51)  yields  Q ^  =  <j>. 

2.  When  pt  =  8,  Equation  (41)  gives  at  =  -<j>.  Then  Equation  (61)  reduces  to  ( h.r)8 , 

~  s 

and  Equation  (51)  yields  Qp  =  8. 

As  a  result 


(64) 


One  motivation  for  introducing  r.h  in  Equation  (51)  was  to  obtain  a  unit-length  Q 


Pt 


Equation  (64)  indicates  that  the  “return”  from  each  individual  SB  interaction  is  co-polarized 
with  the  transmitted  wave,  both  for  horizontal  transmit  (<£)  and  for  vertical  transmit  (8). 


3.6  Structures  of  the  Illuminated  Regions  S12  and  S21 

The  problem  of  determining  the  region  of  integration  Sab  appearing  in  Equation  (52)  must 
be  examined  carefully,  especially  because  it  appears  to  have  been  mishandled  in  the  literature. 

The  region  Sab  is  found  by  tracing  all  the  rays  incident  on  5.4  through  their  reflections  on 
Sa •  Each  point  of  Sb  reached  (i.e.,  illuminated)  by  a  ray  reflected  from  SA  is  part  of  Sab ■  Note 
that  this  section  is  based  exclusively  on  ray  tracing,  i.e.,  on  techniques  from  GO. 

Each  incident  ray  travels  along  direction 


f  —  (  lyi  lz) 


(65) 
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and  goes  through  the  plane  of  SA  at  some  point  xA.  If  xA  is  on  SA,  then  the  ray  is  reflected  along 
direction  [see  Equation  (45)  with  both  sides  divided  by  k)] 

kAB  —  ki  —  kt  —  2(kt-nA)hA 

and,  thus,  goes  through  the  plane  of  SB  at  some  point  xB  given  by 


xB  =  xA  +  aAB  kAB 


(67) 


where  (xAB  is  to  be  determined. 

In  specializing  these  results  to  the  S\S2  and  S2S1  interactions,  it  is  useful  to  temporarily 
restrict  the  analysis  to  the  case 


k  <  1 , 1,  >  0. 

lx 


(68) 


Of  course,  we  are  only  interested  in  the  case  where  lx,ly  >  0.  Afl  these  constraints  are  equivalent 
to 


0  <  <j>  <  7r/4  and  0  <  0  <  7r/2. 


(69) 


These  restrictions  will  be  removed  later. 


3.6.1  S1S2  Interaction 

Here,  one  has  hA  =  y  and  thus  Equation  (66)  can  yield  k  12 
can  be  written  as 


/  0  ^ 

(  Xl  \ 

'  ~lx  ^ 

V2 

= 

0 

+  OL 12 

ly 

V  z2  ) 

l  } 

V  -h  / 

(-lx,ly,-lz).  Equation  (67) 


(70) 
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Eliminating  a  i2  and  solving  for  (2 1,21),  one  gets 


J  xi 

I  zi 


(71) 


By  imposing  the  constraint  that  the  intercept  point  (an,zi)  be  in  Si  (in  order  for  the  incident  ray 
to  be  reflected),  i.e.,  that 


0  <  Xi  <  a  and  —  |  <  z\  <  | 
one  obtains  the  following  constraints  on  (x2,z2) 


(72) 


(73) 


In  addition,  (2:2,22)  must  be  in  S2  (in  order  for  the  reflected  ray  to  be  scattered  by  S2 ),  so  that 


I  0  <  V2  <  a  ^ 

1  “I  <  *2  < 

These  two  sets  of  inequalities  fully  define  S\ 2. 

Careful  analysis  of  these  constraints  reveals  that,  under  the  limitations  of  Equation  (68),  two 
canonically  different  structures  for  region  5j2  must  be  considered.  The  first,  5^,  corresponds  to 
the  case  lz/lx  <  b/a ,  and  the  second,  5 f2,  to  lz/lx  >  b/a  (see  Figure  3).  This  dichotomy  appears 
to  have  been  overlooked  in  the  literature  (a  more  detailed  discussion  is  given  later). 
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3.6.2  S2S1  Interaction 


Here,  one  has  ua  =  x  and  thus  Equation  (66)  yields  k2i  =  (lx,  —ly,  —4)-  Equation  (67)  can 
be  written  as 


f  0  1 

(  l  \ 

lx 

0 

= 

2/2 

+  «21 

—  ly 

(75) 

V  J 

\  22  j 

\  ) 

Eliminating  a2 1  and  solving  for  (y2,  z2),  one  gets 


V2 

Z2 


ly 

4Xl 

Z\  +  jZ-Xl 


(76) 


By  imposing  the  constraint  that  the  intercept  point  (y2,z2)  be  in  S 1  (in  order  for  the  incident  ray 
to  be  reflected),  i.e.,  that 


0<y2  <  a  and  —  f  <  z2  <  f 
one  obtains  the  following  constraints  on 


(77) 


0  <  xi  <  alf 

ly 

<  zi  +  lfx  1  +  |  >  0  (78) 

Zi  +  —  I  ^  0. 

In  addition,  (»i,  sq)  must  be  in  5i  (in  order  for  the  reflected  ray  to  be  scattered  by  Si),  so  that 


0 

b 

2 


< 

< 


Xi 

Z\ 


< 

< 


a 

b 

2 


(79) 


These  two  sets  of  inequalities  fully  define  S2 1- 

Here  also,  two  canonically  different  structures  for  region  S2\  must  be  considered.  The  coun¬ 
terpart  of  Figure  3  is  Figure  4. 
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3.7  Simple  Expressions  for  1^2  an(i  pS2S  1 

Having  determined  the  structures  of  the  regions  5i2  and  £21,  one  can  now  evaluate  pSiS2 
and  P^^Si  using  Equation  (52). 

3.7.1  SiS2  Interaction 

With  fiA  =  y,fiB  =  x ,  and  r'  =  (0,^,z),  Equation  (52)  can  be  written  as 


=  lxj  e™l*zdS12. 


pSiS2  —  1  I  „2iklzZ 


IS  12 


(80) 


When  lz/lx  <  b/a ,  the  region  5^2  (see  Figure  3)  must  be  used  and  Equation  (80)  becomes 

(81) 


p5i£2 


ly 

rt 

r  b  l  1 

r  2--ry  , 

/  .  lv  e2ikl*zdz 

Jo 

~  2 

dy. 


When  lzllx  >  b/a,  the  region  S12  (see  Figure  3)  must  be  used  and  Equation  (80)  becomes 

dy.  (82) 


ft 

f2  e2iklzzdz 

Jo 

J-2 

3.7.2  S2Si  Interaction 

With  fiA  =  x,fiB  =  y ,  and  r'  =  (x,0,z),  Equation  (52)  can  be  written  as 


pS2Sl  =  ,’Js 


e2lkl*zdS21. 


(83) 
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For  lzjlx  <  b/a,  one  finds  (see  Figure  4) 


PS^=l,£  ftf‘‘ 


dx. 


(84) 
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For  lz/lx  >b/a,  one  finds  (see  Figure  4) 


Ps2s1 


.  f  bY  f  2 ' 

”  yL  z  U 


,b_l_z_x 

2  lx  2iklzz 

b  e 


*zdz 


dx. 


(85) 


It  is  useful  to  observe  that  the  change  of  variable  y  =  y'(ly/lx )  in  Equations  (81)  and  (82) 
yields  Equations  (84)  and  (85),  respectively.  In  other  words,  in  all  cases, 


pS\  S2  _  p>?2  S\ 


(86) 


This  equality  can  also  be  proved  graphically.  It  is  clear  from  Figures  3  and  4  that  S\2  is  simply  52i 
“compressed  horizontally”  by  a  factor  ly/lx  (note  the  common  role  of  z0  in  the  figures).  Therefore, 
the  ratio  of  the  integrals  in  Equations  (80)  and  (83)  is 


!S„e*t,“dSn  l, 
fs21  e2tk^zZdS2i  b 


(87) 


Equation  (86)  follows  immediately.  Note  that  this  equality  is  partially  a  result  of  the  inclusion  of 
r.fiB  in  Equation  (52). 

Straightforward  integration  of  Equations  (84)  and  (85)  yields 


pS1S2  =  pSiS!  = 


dbi 


’ly 

2ikblz 

ably 

2ikdxlz 


,ik(b— dx)l. 


sine  k'blz  —  e 


'z  sine  k'dxlz  —  e 

~zkbl z 


—ikbl 


■] 

if  k 


(88) 


where 


dT  =  a- 


(89) 


As  one  would  expect,  the  results  of  Equation  (88)  become  identical  when  lz/lx  =  b/a,  in  which  case 
dx  =  b. 
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3.8  Simple  Expressions  for  Qp^2  and  Qp^1 

_  c  c 

One  begins  by  rewriting  the  numerator  B  of  the  Equation  (53)  for  B  as 

B  =  [I-  rf](nB  X  [at  -  2(at.nA)hA\)  (90) 

=  [I  -  fr](hB  x  at)  -  2(at.hA)[I  -  rr](hB  x  hA)  (91) 

=  ( hB.(at  x6)~  2 (at.hA)[0.(hB  x  hA)})0  +  hB.(at  x  4>)4>  (92) 


qSaSe  _ 

<i> 


=  -2- 


where  one  has  expanded  the  first  and  second  [/-ff](...)  expressions  in  Equation  (91)  by  expressions 
similar  to  Equations  (63)  and  (62),  respectively,  and  used  the  fact  that  hB  X  hA  is  necessarily 
orthogonal  to  4>. 

1.  When  pt  =  Equation  (41)  gives  at  =  0.  Using  this  value  for  at  and  noting  that 
f  =  6  x  4>,  Equation  (92)  can  be  reduced,  and  Equation  (53)  yields 

(93) 

r.nB 

To  further  reduce  the  fraction  above,  proceed  as  follows.  First,  because  nB  X  hA  is 
parallel  to  i,  one  can  express  this  product  solely  in  terms  of  its  components  along  0 
and  f,  i.e., 

nBxhA  =  0[d.(nB  x  n^)]  +  r[f.(hB  x  hA)].  (94) 

Second,  transform  each  side,  say  .s,  of  this  expression  by  the  operation  s.hA,  and  use 
the  fact  that  ( nB  X  hA).hA  =  0.  One  finds 

(§.hA)[§.(hB  X  hA)\  +  ( f.hA)[r.(hB  x  hA )]  =  0.  (95) 

Finally,  extracting  from  this  expression  the  numerator  appearing  in  Equation  (93), 
one  gets 


qSASB  =  2{r.nA)[f.(hB  X  hA)]£  +  ^ 


(96) 

(j)  r.nB 

For  the  S1S2  interaction,  one  has  nA  =  y  and  hB  =  x,  so  that  Equation  (96)  reduces 
to 


QS'S*  =2lflJ  +  j>. 

<p 

For  the  S2S1  interaction,  one  has  nA  =  x  and  fiB  =  y,  so  that 
qS2si  =  _2  l±lz§  +  l 

q>  ly 


(97) 


(98) 
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(99) 


2.  When  pt  =  6,  Equation  (41)  gives  at  =  -<j>,  and  then  Equation  (53)  yields 

qSaSb  =(1  +  J^A)[0-(nB  X  nA)]\  ~ 

^6  \  f.fiB  ) 

To  further  reduce  the  fraction  above,  proceed  as  follows.  First,  because  <j>  is  orthog¬ 
onal  to  z,  one  can  express  <p  solely  in  terms  of  its  components  along  nA  and  hB, 
i.e., 

cf>  =  fiA^-fiA)  +  nB(4>-nB)-  (10°) 

Second,  transform  each  side,  say  s,  of  this  expression  by  the  operation  0.(hB  X  s). 

One  finds 

0.{tib  x  <j>)  =  (4>.hA)[0.(nB  x  ua)\  (101) 

or,  using  the  previously-used  cyclic  property  of  the  mixed  product, 

{^.nA)[0.{hB  X  hA)]  =  -ns-r.  (102) 

Finally,  extracting  from  this  expression  the  numerator  appearing  in  Equation  (99), 
one  gets 

qSaSb  _  _q  (103) 

6 

and,  therefore, 

QSJS 2  =  QSJS 1  =  -6.  (104) 

6  6 

One  motivation  for  introducing  r.fiB  in  Equation  (53)  was  to  obtain  unit-length 

QSJS2  and  Q^2^1  (the  co-polarized  components). 

0  0 

For  horizontal  transmit  (<£),  Equations  (97)  and  (98)  indicate  that  the  “return”  from  each 
individual  DB  interaction  has  components  that  are  both  co-polarized  and  cross-polarized  with  the 
transmitted  wave.  However,  when  the  sensor  is  in  the  plane  (x,y),  where  lz  =  0,  the  cross-polarized 
component  vanishes.  Further  conclusions  regarding  this  component  will  be  given  when  the  joint 
return  for  the  SiS2  and  5’25'i  interactions  are  considered. 

In  the  case  of  vertical  transmit  ((?),  Equation  (104)  indicates  that  the  “return”  from  each 
individual  DB  interaction  is  totally  co-polarized  with  the  transmitted  wave. 
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3.9  Final  Expression  for  the  Backscattered  Fields 

Some  of  the  results  of  Sections  3.6  and  3.7  were  derived  under  the  limitations  of  Equation  (68), 
i.e.,  I  <  lx  and  4  >  0.  Simple  symmetry  arguments  (or  brute-force  recalculation)  show  that  all 
results  can  be  generalized  to  arbitrary  values  of  lx,  ly,  and  lz  (in  the  allowed  ranges)  by  defining 


I  =  max(4,4) 
m  -  min(4,4) 

TO  =  |4| 


(105) 


and  making  the  substitutions 


m 


(106) 


throughout  the  preceeding  results. 

In  particular,  the  important  dichotomy  related  to  the  sign  of 


4  _  6 

4  a 


can  be  generalized  by  considering  the  sign  of 


n 


14 


(107) 


(108) 


l  a  max(4,  ly)  a 
Because  of  the  observed  equalities  in  Equations  (64)  and  (86),  one  can  rewrite  Equation  (49) 


as 


EpSr^  2tt  r 


or,  in  terms  of  SB  and  DB  contributions, 


ikE0e  ikT  (pSBpSB  pDBnDB) 

Ept(r)  - - 2^7" QPt  +  ^Pt  } 


=  JkE°e  i  +  p52)^1  +  P^HQlf2  +  Qpf1)} 


(109) 


(110) 
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where 


pSB 

qsb 

Vpt 

pDB 

q?b 

y Pt 


pS i  _|_  pS2 

QSJ 
* Pt 

2pSi$2 


=  0.5  (Q 


r>SlS2  ,  /y^l 


Pt 


+  Q-‘ 

T  v Pt 


(111) 


The  reason  for  the  seemingly  arbitrary  factors  2  and  0.5  above  is  explained  below. 

The  values  of  the  four  fundamental  quantities  in  Equation  (111)  can  be  expressed  by  making 
the  substitutions  of  Equation  (106)  in  Equations  (57),  (60),  (64),  (88),  (97),  and  (98). 

One  finds 


PSB  =  ab  sine  k'bn(m  eikal  sine  k'al  +  l  elkam  sine  k' am ) 


(112) 


=  Pt  = 


4>  if  pt  -  4>  (H  pol) 
6  if  fit  =  0  (V  pol) 


fp^(etk(b-din)n  sine  k'dinn  —  e  lkbn )  if  din  <  b 
-^_(sinc  k'bn  -  e~ikbn )  if  dtn  >  b 

ikdlnn  v 

j  -Clmnd  +  4>  if  pt  =  <t>  (H  pol) 

I  -6  if  Pt  =  0  (V  pol) 


with 


and 


din  —  ®  ^ 


^■Imn  — 


l 2  —  m2 
Im  71 


(113) 


(114) 


(115) 


(116) 


(117) 
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One  motivation  for  introducing  the  extra  factors  2  and  0.5  in  Equation  (111)  was  to  obtain  unit- 
length  vectors  for  the  co-polarized  components  in  Equation  (115). 

3.10  Backscattered  Field  for  kt  Orthogonal  to  Crease 

In  the  special  case  where  the  transmitted  wave  travels  in  a  direction  orthogonal  to  the  crease 
of  the  dihedral  (i.e.,  when  n  =  0),  Equations  (112)  to  (117)  can  be  reduced  considerably.  All 
simplifications  are  straightforward,  except,  perhaps,  in  the  case  of  Equation  (114).  Because  din  =  0, 
the  first  expression  for  PSB  must  be  used.  Also,  in  the  limit  as  n  0,  dinn  approaches  0  as  n2, 
whereas  bn  approaches  0  as  n.  Therefore, 

PDB(n  =  0)  =  lim  ^( eikbn  -  e~ikbn ) 
n-o  ikbny  ’ 

=  2abm  lim  sine  k'bn 

n— >0 

=  2  abm.  (118) 


Thus,  for  n  —  0, 

PSB(n  =  0)  =  ab  (m  eikal  sine  k'al  +  l  eikam  sine  k'am ) 


(119) 


QSB(n  =  0  )  =  pt  = 


<i>  if  pt  =  4>  (H  pol) 
9  \i  pt  —  9  (V  pol) 


(120) 


PDB{n  =  0)  =  2  abm 


(121) 


QDB(n  =  0)  = 


^  if  pt  =  4>  (H  pol) 
-6  iipt  =  e  (V  pol) 


(122) 


with,  of  course, 

l  —  \/l  —  rn2  or  m  =  \/l  —  l2. 


(123) 
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4.  DISCUSSION  OF  RELATION  BETWEEN  n/l  AND  b/a 


This  section,  explores  the  origin  of  the  dichotomy  related  to  the  sign  of  the  quantity  n/l-  b/a 
[see  Equation  (108)],  and  analyzes  the  mapping  of  this  sign  on  a  3-D  viewing  sphere. 

4.1  Fundamental  Observations 

Consider  the  dihedral  positioned  in  the  (x',y',z')  coordinate  system  as  indicated  in  Figure  5 
(note  that  the  origin  is  at  the  bottom  end  E  of  the  dihedral’s  crease).  To  simplify  the  discussion, 
consider  the  special  form  lz/lx  -  b/a  of  the  quantity  n/l  -  b/a.  As  previously  discussed,  this  special 
case  arises  when  ly  <  lx  (i.e.,  0  <  <p  <  tt/4)  and  lz  >  0  (i.e.,  0  <  6  <  x/2).  The  conclusions  reached 
for  this  special  case  are  easily  extended  to  the  general  case  via  symmetry. 

Consider  the  plane  AGE  that  contains  the  diagonal  AE  of  S\  and  is  perpendicular  to  Si.  It  is 
simple  to  show  that  any  transmitted  plane  wave  with  wavevector  kt  —  (lx,ly,lz)  in  this  plane  (and 
with  ly  <  lx  and  lz  >  0)  satisfies  the  equality  of  interest,  i.e.,  lz/lx  =  b/a.  Indeed,  from  Figure  5, 
El  =  b tan#  and  El  —  a/  cos <f>,  so  that 


cot  6  b  lz  b 

T  —  ~  >  7~  ~  ~ 

cos  <pa  lx  a 


(124) 


where  Equation  (55)  was  used. 

If  one  (a)  selects  a  particular  wavevector  in  the  previously  considered  plane  AGE ,  say  with 
direction  GE  (the  corresponding  kt  is  shown  in  Figure  5),  (b)  considers  the  plane  GCE  containing 
GE  and  perpendicular  to  S2,  and  (c)  computes  the  y'  coordinate  of  the  intersection  C  of  that  plane 
with  GH  or,  equivalently,  the  y'  coordinate  of  F,  one  finds  y'F  =  El  sin  <p  =  a  tan  cp  —  aly/lx.  Note 
that  the  line  containing  CF  is  also  the  line  appearing  in  S2  in  Figure  3. 

This  discussion  shows  that  any  wavevector  kt  —  (lx,ly,lz)  satisfying  lz/lx  =  b/a  (and  ly  <  lx, 
lz  >  0)  is  defined  by  the  intersection  of  two  planes:  (a)  the  plane  perpendicular  to  S\  and  containing 
the  diagonal  AE  of  Si,  and  (b)  the  plane  perpendicular  to  S2  and  containing  the  diagonal  GE  of 
the  rectangle  BCFE  uniquely  defined  by  aly/lx.  Similar  conclusions  can  be  reached  for  conditions 
other  than  ly  <  lx  and  lz>  0. 

The  previous  reasoning  can  be  carried  out  elegantly  via  spherical  trigonometry.  Using  the 
sphere  shown  in  Figure  5  (with  center  E  and  radius  a),  consider  some  viewing  direction  GE.  The 
intersection  0  of  this  fine  with  the  sphere,  together  with  the  three  coordinate  axes  for  x',  y' ,  and 
z',  define  the  arcs  of  great  circle  DOM,  JOK,  and  LON,  respectively  (e.g.,  the  plane  through  0 
and  the  x'  axis  defines  DOM). 
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By  applying  one  of  Napier’s  rules  [21,  pp.  271-272]  to  the  right  spherical  triangle  ONJ,  one 

finds 


cos  <j>  —  cot  8  cot  ay  (125) 

which  implies  that  lz/lx  =  b/a  can  be  satisfied  only  if  tana,,  =  b/a,  i.e.,  only  if  the  plane  JKE 
contains  the  diagonal  AE  of  S\. 

By  applying  a  similar  formula  to  OND,  one  finds 


sin  <j>  =  cot  6  cot  ax 


(126) 


which,  given  the  previous  constraint  on  ay,  gives  tana^  =  ( b/a)(lx/ly ),  and  thus  EF  —  aly/lx. 

4.2  Partitioning  of  the  Viewing  Sphere 

Another  advantage  of  the  analysis  on  the  sphere  is  that  it  provides  an  instant  display  of  the 
partitioning  of  the  viewing  space  by  the  boundary  case  lz/lx  =  b/a.  Indeed,  the  boundary  on  the 
sphere  of  Figure  5  is  the  arc  KOJ ,  restricted  to  0  <  <j>  <  7r/4  because  of  the  current  limitation  to 
ly  <  lx  and  lz  >  0. 

A  simple  symmetry  argument  readily  extends  the  conclusions  to  the  general  boundary  case 
n/l  =  b/a.  Figures  6(a)  and  6(b)  show  the  regions  on  the  sphere  that  correspond  to  n/l  <  b/a  and 
n/l  >  b/a  for  a  long  and  a  short  dihedral,  respectively  (only  for  lz  >  0). 

4.3  Physical  Interpretation 

When  lz/lx  -  b/a ,  it  is  clear  from  Figures  3  and  4  that  F'  merges  with  F  and  D’  merges  with 
D.  The  resulting  diagonals  BF  and  BD  are  shown  in  Figures  7(a)  and  7(b). 

When  lz/lx  =  b/a,  Equation  (71)  reveals  that  the  ray  with  direction  kt  that  is  reflected  at  A 
by  Si  arrives  at  F  on  S2.  It  is  thus  clear  that  the  sheet  of  rays  that  is  reflected  by  the  edge  AB 
exactly  illuminates  the  diagonal  BF,  i.e.,  all  rays  incident  along  AB  hit  S2  (along  BF).  This  is 
illustrated  in  Figure  7(a). 

As  previously  observed,  the  direction  kt  must  be  parallel  to  both  the  plane  AEF  and  the 
plane  CED.  Furthermore,  Equation  (76)  reveals  that  the  ray  with  direction  kt  that  is  reflected  at 
C  by  S2  arrives  at  D  on  Si-  Clearly,  the  sheet  of  rays  reflected  by  the  partial  edge  BC  exactly 
illuminates  the  diagonal  BD.  This  is  illustrated  in  Figure  7(b). 
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4.4  Interpretation  of  Parameter  dx 


A  simple  interpretation  of  the  parameter  dx  defined  by  Equation  (89)  can  be  obtained  from 
Figure  8  (which  has  many  elements  in  common  with  Figure  5).  As  shown  in  Section  4.1,  the 
arc  JOK  (where  K  is  on  the  diagonal  AE)  corresponds  to  all  incidence  directions  kt  related  by 
lz/lx  =  b/a.  This  observation  and  the  definition  in  Equation  (89),  i.e.,  lz/lx  =  dx/a,  immediately 
suggest  a  graphical  construction  for  dx.  By  first  drawing  the  arc  JO'K'  through  the  point  O' 
corresponding  to  the  direction  kt  of  interest,  and  then  finding  the  intersection  A'  of  the  lines  EK' 
and  AD,  one  finds  that  dx  is  simply  the  length  of  the  segment  DA'.  Clearly,  an  important  value 
for  dx  [and  its  equivalent  din  given  by  Equation  (116)]  is  6,  which  corresponds  to  A'  =  A.  Indeed, 
depending  upon  the  relative  values  of  dx  (or  din)  and  6,  different  expressions  must  be  used  for  the 
DB  contribution  given  by  Equation  (88)  [or  Equation  (114)]. 
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Figure  5.  Geometry  for  studying  sign  of  lz/lx  —  b/ a. 
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e  6.  Partitioning  of  the  viewing  sphere  according  to  the  sign  ofn/l  —  h/a  for  ( 
lihedral  (large  b/a)  and  (b)  a  short  dihedral  (small  h/a). 
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5.  POLARIZATION  SCATTERING  MATRIX  OF  DIHEDRAL 


As  indicated  in  Section  2.4,  one  should  first  derive  the  complex  radar  cross  sections 
given  by  Equation  (15).  Then  one  can  obtain  the  elements  Sp^  of  the  scattering  matrix  through 
Equation  (16).  The  corresponding  matrices  Jo  and  S  can  also  be  related  through  Equation  (20). 

To  continue  keeping  single-bounce  (SB)  and  double-bounce  (DB)  contributions  distinct,  Equa¬ 
tion  (110)  is  rewritten  as 


Ep,(f)  =  £ff  (f)  + 

with 

ES-B(r)  -  J^rUpSBQSB 
~  2  xr  QP* 

)  _  JkE0e-ikT  dbqDB 

LPt  ^  _  2 xr  ' 

Then,  Equation  (15)  can  be  rewritten  as 


(127) 


(128) 

(129) 


(130) 


with 


lim  2jrr 

r— t-oo 


hm  2jirr 

r— »■  oo 


ZSpf(r).pr 

E« 


Eo 


(131) 

(132) 


5.1  Single-Bounce  Complex  Radar  Cross-Section 

The  general  expression  for  the  SB  C-RCS  is  obtained  from  Equations  (131),  (128),  (112), 
and  (113),  i.e., 


SB  - 

Qpt  -Pr 


K™  Pt.Pr 


(133) 
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where 


K, 


SB 


Therefore 


=  KSB 


(134) 


(135) 


5.2  Double-Bounce  Complex  Radar  Cross-Section 

The  general  expression  for  the  DB  C-RCS  is  obtained  from  Equations  (132),  (129),  (114), 
and  (115),  i.e., 


where 


DB 


r\DB  „ 

Qpt  •Pi' 


=  K?B 


Iran  O-Pr  +  <t>-Pr  If  Pt  =  4>  (H  pol) 
-6-Pr  if  Pt  =  0  (V  pol) 


DB 


Therefore 


1 


€/ran 


0 

-1 


(136) 


(137) 


(138) 


where  [see  Equation  (117)] 


£lmn 


l2  -  m2 
Im  n 


(139) 
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For  an  arbitrary  triplet  of  consistent  values  /,  m,  and  n  (e.g.,  requiring  that  l  >  m),  eimn 
is  generally  nonzero  and  thus  the  matrix  \J oDB  is  not  symmetric.  When  l  >  m,  one  has  e/mn  « 
(/ /m)n,  and  this  value  can  be  large  if  n  is  not  too  small.  However,  these  conditions  correspond  to 
large  deviations  of  the  wavevector  kt  from  the  principal  axis  of  symmetry  (characterized  by  l  =  m 
and  n  =  0)  and  thus  to  viewing  directions  where  the  PO  approximation  would  not  be  expected 
to  produce  accurate  results.  For  example,  it  is  reported  in  Kouyoumjian  [7,  p.  871]  that  the  PO 
prediction  for  scattering  from  a  circular  plate  with  radius  a  is  satisfactory  only  within  20  deg  of 
normal  incidence  for  ka  >  8.5. 

However,  for  orientations  kt  parallel  either  to  the  plane  l  =  m  or  to  the  plane  n  =  0,  Qmn  is 
precisely  zero.  Furthermore,  eimn  is  very  small  in  the  vicinity  of  the  principal  axis  of  symmetry, 
where  PO  is  expected  to  produce  accurate  results.  As  a  result,  one  sets  e;mn  to  zero  in  subsequent 
discussions,  thereby  making  V<JDB  symmetric,  i.e., 

^  j  .  (140) 


The  reason  for  the  presence  of  eimn  in  Equation  (138),  and  thus  the  lack  of  symmetry  of 
yJoDB ,  is  the  failure  of  PO  to  obey  the  reciprocity  theorem  [7,  8],  with  the  result  that  bistatic 
(double-bounce  in  this  case)  cross  sections  are  generally  erroneous. 

5.3  Relative  Importance  of  Single-  and  Double-Bounce  Complex  Radar  Cross  Sec¬ 
tions 

To  compare  the  magnitude  of  the  SB  and  DB  C-RCS  near  the  main  axis  of  symmetry  of 
the  dihedral,  we  examine  the  ratio  of  the  factors  KBB  and  KBB  [appearing  in  Equations  (135) 
and  (138)]  in  the  case  where  l  =  m  =  V2/2  and  n  =  0.  Using  Equations  (119)  and  (121),  one  finds 


P  = 


Ksb 


KDB 

cry  sym 


pSB 

1  sym 


pDB 

M  sym 


•  y  ^ 

smc  k  a~ 


(141) 


Observe  that  p  decreases  towards  0  as  ka  increases.  Furthermore,  assuming  ka  >  8  (the  requirement 
used  in  Kouyoumjian  [7,  p.  871]  for  PO  to  produce  accurate  results),  one  has 

=  —  «  0.212.  (142) 

1.5  TT 


P  < 


sin  1.57T 


1.57T 
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Therefore,  the  SB  C-RCS  can  generally  be  assumed  to  be  at  least  13.5  dBs  below  the  DB  C-RCS. 
Even  though  this  result  was  derived  on  the  main  axis  of  symmetry,  there  is  a  tendency  for  the  DB 
C-RCS  to  dominate  whenever  n  is  small.  Keeping  in  mind  that  the  results  obtained  through  the 
PO  approximation  are  valid  only  in  the  vicinity  of  the  main  symmetry  axis,  the  contribution  of  the 
SB  C-RCS  can  generally  be  ignored,  and  the  C-RCS  matrix  of  the  dihedral  is  often  taken  as 

^~V^DB  =  KDB  ^  °^j.  (143) 

The  observation  that  the  DB  interaction  is  generally  more  important  than  the  SB  interaction  was 
also  made  in  Corona  et  al.  [10]. 

5.4  Double-Bounce  Complex  Radar  Cross  Section  on  the  Symmetry  Axis 

Setting  l  =  m  =  V2/2  and  n  =  0,  and  using  Equations  (137)  and  (121),  one  finds 


K, 


DB 

cr,sym 


=  — ^g6v/2. 


(144) 


The  co-polarized  DB  radar  cross-sections  (RCSs)  are  thus  found  to  be 

2  k2a?b2  8?r  a2b2 


(145) 


° 4>4>  ° 99  7r  A2 

These  results  agree  with  the  RCS  values  published,  e.g.  in  Knott  [11,  Table  6-1,  p.  178]. 

5.5  C-RCS  Matrix  in  Circular  Basis 

By  applying  the  equivalent  of  Equation  (31)  for  C-RCS  matrices,  one  can  derive  from  the 
dihedral  linear-basis  C-RCS  matrix  [see  Equation  (143)] 


1  0 

0  -1 


(146) 
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the  dihedral  circular-basis  C-RCS  matrix 


y/°(LK)  —  K?B 


(147) 


This  last  expression  [and  not  Equation  (146)]  makes  the  dihedral  an  even-bounce  (circular  polar¬ 
ization)  scatterer. 


45 


APPENDIX  A 

TRANSFORMATION  OF  SCATTERING  MATRIX  UNDER  CHANGE  OF 

POLARIZATION  BASIS 

The  rules  governing  the  transformation  of  the  scattering  matrix  when  a  change  of  polarization 
basis  is  performed  are  well-known.  They  are  included  here  for  completeness.  This  report’s  approach 
is  a  mix  between  those  of  Boerner  et  al.  [22,  23]  and  Riegger  [24] . 

A.l  Transformation  of  Basis  Vectors 

Consider  the  transformation  between  the  polarization  bases  ( AB )  and  ( A'B ')  respectively 
characterized  by  their  (unit-length)  basis  vectors  hA,hB  and  hAi,hBi,  and  assume  that  these  vectors 
are  related  by  Equation  (24),  i.e., 

h(A'B')  —  Th{AB)  (A-l) 

or 

(  ^A‘  \  _  (  Ta'a  Ta'b  \  /  hA  \ 

\  hB'  J  \  TB’A  TB:B  )  \hB  ) 

As  one  shall  see  shortly,  T  is  necessarily  a  unitary  matrix. 

A. 2  Transformation  of  Fields 

Any  field  E  can  be  expressed  in  each  of  the  (AB)  and  ( A'B ')  bases  as 

E  =  EAhA  +  EBhB 
E  =  EA'h.A<  +  EB'hB>. 

By  using  Equation  (A.2),  Equation  (A.4)  becomes 

E  -  EA'(TA'AhA +  TA'BhB)  +  EB'(TB'AhA +  TB'BhB) 

=  (TAiAEAi  +  TBiAEBi)hA  + (TA'BEA>  +  TBiBEB’)hB.  (A. 5) 


(A.3) 

(A.4) 
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Direct  comparison  of  Equation  (A. 3)  and  (A. 5)  shows  that 


I  EA  \  =  /  TA'a  Tb,a  W  Ea<  \  (A .6) 

\  Eb  J  \  TA'b  TB'b  /  \  Eb>  ) 

or,  more  compactly, 

E(ab)  -  UE(a'B')  (A'7) 

where,  in  general,  E^xy)  denotes  the  two-element  vector  containing  the  components  of  E  in  the 
(AY)  basis,  and 

U  =  Tt.  (a-8) 

The  fact  that  the  total  power  in  the  (plane)  wave  corresponding  to  E  or,  equivalently,  the  magnitude 
of  E ,  must  remain  constant  under  a  change  of  basis  leads  to  the  requirement  that  U ,  and  thus  T, 


be  unitary  matrices,  i.e., 

U~l^U*T  and  |  det(?7)|  =  1  (A-9) 

T~1=T*t  and  |  det(T)|  =  1.  (A-10) 

If  ( AB )  is  the  linear  basis  (HV):  and  (A'B')  an  elliptical  basis  with  arbitrary  orientation 
and  arbitrary  ellipticity  r,  one  can  show  that  [24] 

E(hv)  =  U(i?,T)E(A'B')  (A-11) 

with 

U(iP,t)  —  R(iP)H(t)  (A.12) 

where  R(ib)  is  the  rotation  matrix  . 


W)  = 


cos  ip  —  sin  ip 
sin  ‘ip  cos  xp 


(A.13) 
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and  H (r)  the  ellipticity  matrix 


H(t)  = 


cos  r  j  sin  r 
j  sin  r  cos  r 


(A- 14) 


It  is  useful  to  note  that  R(ip)  and  H(r)  are  unitary,  with  the  consequence  that  their  product,  i.e., 
t/(V>,r),  is  also  unitary  (as  previously  announced). 

A. 3  Transformation  of  Scattering  Matrices 

The  starting  point  is  the  Equation  (9)  defining  the  scattering  matrix  S  in  the  linear  basis 
(HV),  i.e., 


E\hv)  =  a(k,'r)S(HV)Et(HV)-  (A.15) 

Now,  as  strange  as  it  may  seem,  one  will  perform  distinct  changes  of  basis  on  the  received  and 
transmitted  fields!  ET  is  expressed  in  a  basis  (A'rB'r)  characterized  by  ip  and  rr,  and  Et  is  expressed 
in  a  basis  ( A'tB[ )  characterized  by  w  and  rt.  Therefore,  using  Equation  (A.ll),  Equation  (A.15) 
becomes 

U(^,TT)Er(A>rBi^  =  a(k,r)S(HV)U{‘ip,Tt)Et(AitB't)  (A.16) 

or 

E>\A’TB’r)  =  akjU-l(i>,Tr)S(HV)U{'<J),Tt)Et(A'tB'ty  (A.17) 

By  generalizing  Equation  (A.15)  to  the  case  of  distinct  receive  and  transmit  bases,  one  can 

write 


^(A'rB'r)  =  a(k^S{A!rB'rA[B't)E\A'tB'ty 


Comparing  Equations  (A.17)  and  (A. 18),  one  gets 


S 


(A'TB'r,A'tB't) 


(A.18) 


(A. 19) 
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When  Tr  —  rt  =  r,  Equation  (A. 19)  describes  a  true  change  of  basis ,  i.e., 


S(A'B')  =  u  ls(HV)U  =  U*TS(HV)U  (A.20) 

where  {MB')  =  ( A'tB[ )  =  ( A'rB'r )  and  U  =  U(i/>,t).  However,  this  is  not  the  formula  to  use  for 
transforming  5. 

The  correct  transformation  is  obtained  by  using  rr  =  -rt  =  -r,  in  which  case  Equation  (A.19) 
describes  a  pseudo  change  of  basis ,  i.e., 


5 


U  l(ip,-T)S(HV)U{ip,T). 


(A.21) 


It  is  easy  to  see  that  U  ~T)  —  so  that 

S(A'B')  =  UtS(HV)U  (A.22) 

where  it  is  understood  that  U  =  U(if,T)  corresponds  to  the  (AtBt)  basis  and  that  the  subscript 
(A'B')  corresponds  to  the  basis  the  transmitted  field  is  transformed  into.  The  transformation 
applied  to  the  received  field  need  not  be  known  in  order  to  use  Equation  (A.22),  but  it  is  easy  to 

see  that 

J7(^,-r)  =  {UT)~l  =  U*.  (A-23) 

Even  though  Equation  (A.22)  is  the  classical  formula  for  transforming  a  scattering  matrix, 
it  is,  once  again,  important  to  understand  that  this  formula  does  not  describe  a  true  change  of 
basis  [24,  p.  21]. 

Equation  (A.22)  can  easily  be  generalized  to  handle  transformations  between  arbitrary  bases, 
say  ( AB )  and  (A'B').  Using  Equation  (A.22),  one  has 


S(AiBi)  =  UlS(HV)Ul  (A.24) 

S(A2B2)  =  UjS{HV)U2. 
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It  follows  that 


(A.25) 

and  thus 

S{A,B,1  =  (Ur'VifS^B^Ur'Ui).  (A.26) 

Furthermore,  Equation  (A. 7)  gives 


{E(HV)  - 

E(HV)  =  ^2-E(A2fi2) 

so  that 

=  UilU2E{MB2y 

Therefore,  it  is  clear  from  Equations  (A.26)  and  (A. 28)  that  if 


(A.27) 


(A.28) 


E(ab)  =  UE(A'B’)  (A.29) 

then 

%'B')  =  UtS(ab)U  (A.30) 

which  is  the  desired  generalization  of  Equation  (A.22).  Furthermore,  note  that  Equation  (A.29)  is 
identical  to  Equation  (A.7)  so  that  U  =  TT  [see  Equation  (A.8)]  where  T  is  the  matrix  relating  the 
basis  vectors  of  the  bases  ( AB )  and  ( A'B '). 
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It  should  be  observed  that  Equation  (A.30)  guarantees  the  conservation  of  the  voltage  (see  [23, 

p.  211]) 


~~T 

V  =  E(AB)S(ab)E{ab) 


(A.31) 


under  a  pseudo  change  of  basis.  Indeed,  using  Equations  (A.29)  and  (A.30),  Equation  (A.31) 
becomes 


V  =  eJAIBi)Ut  S(ab)U  E(a'B') 

_  m  — 

=  E^b^S^b^E^a'B'Y 


(A.32) 
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APPENDIX  B 

ALTERNATE  FORMULAS  FOR  THE  TOTAL  BACKSCATTERED  FIELD 


Here,  one  uses  Equations  (110),  (112),  (113),  (114),  and  (115)  to  derive  a  set  of  new  expressions 
for  the  total  backscattered  field  E^.  These  expressions  are  written  in  a  form  that  is  as  close  as 
possible  to  results  derived  by  Corona  et  al.  [10].  The  algebraic  manipulations  required  to  achieve 
this  goal  have  the  effect  of  mixing  the  SB  and  DB  contributions  (which  may  not  always  be  desirable). 

When  n/l  <  bja ,  one  has 


ikE0e~ikr  , 

— - ab 

2irr 

[sine  k'bn(m  etkal  sine  k'al  +  l  elkam  sine  k'am  ■ 


m 

ikbn 


■eikbn(l  ■ 


e  tkdlnn  sjnc  k'dinn)]4> 


ikE0e~ikT  l2  —  m2 

j - aj) - - - n 

2xr  l 

[2  sine  k'bn  -  7^eife6n(l  -  e~ikd^n  sine  k' 

IKOTl 


(B.l) 


E^if)  = 

ikE0e~ikr  , 

— - ab 

2i rr 

[sine  k'bn(m  elkal  sine  k'al  +  l  elkam  sine  k’am  -  2m) 

+_H!_eikbns1  _  e-ikdinn  si nck'dinn)\e. 
ikbn 


(B.2) 
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When  n/l  >  b/a,  one  has 


)  = 


ikEoe' 


-ikr 


2xr 


-ab 


[sine  k'bn{m  eikal  sine  k'al  +  l  eikam  sine  k'am  +  ) 


lhdlnTl 


e~lkbn](f> 


ikE0e~ikr  l2  —  m2 
-) - f - ab - : - ft 


27rr 


l 


[ — - — (sin ck'bn  —  e  lfc6n)]0 
ikdinii 


(B.3) 


Eg(f)  = 


ikEoe 


-ikr 


2i:r 


-ab 


[sine  k'bn(m  eikal  sine  k'al  +  l  eikam  sine  k’am  - 


I  ra  jkbn^p 
zkdl’n  Tl 


(B.4) 


The  Equations  (B.l)  and  (B.2)  become  very  similar,  but  not  identical,  to  the  Equations  (22) 
and  (21)  of  Corona  et  al.  [10]  when  /,  m,  and  n  are  expressed  in  terms  of  <j>  and  6  through  Equa¬ 
tions  (105)  and  (55).  The  differences  are  limited  to  the  signs  of  some  of  the  terms  in  the  equations 
and  the  signs  of  the  exponents  in  some  of  the  complex  exponentials.  Of  course,  these  sign  differences 
are  important.  The  author  of  this  report  has  not  been  able  to  be  explain  these  differences. 

The  dichotomy  related  to  the  sign  of  n/l  -  b/a  is  overlooked  in  Corona  et  al.  [10].  Ret¬ 
rospectively,  the  analysis  in  [10]  corresponds  to  the  case  n/l  <  b/a.  Thus,  the  equivalents  of 
Equations  (B.3)  and  (B.4)  do  not  appear  in  [13]. 
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